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Abstract 

We study the discovery prospects of a heavy neutral scalar arising from a U(1)b~l extension 
of the Standard Model (SM) during the Large Hadron Collider’s high luminosity runs (HL- 
LHC). This heavy neutral scalar mixes with the SM Higgs boson through a Higgs portal and 
interacts with the SM particles with an interaction strength proportional to the sine of the 
mixing angle. The mixing between the two Higgs bosons is constrained by direct and indirect 
measurements. We choose an experimentally viable mixing angle and explore in detail the ZZ 
and WW decay modes of the heavy Higgs boson. For the ZZ case, we focus on the cleanest 
M and 2l2j final states and find that a heavy Higgs boson of mass smaller than 500 GeV can 
be discovered at the HL-LHC. For the WW decay mode, we analyze the ljj$T signature. We 
implement novel background reduction techniques in order to tackle the huge background by 
performing both cut-based and multivariate analyses. However, large backgrounds render this 
channel challenging. We briefly discuss the discovery prospects of the heavy Z'-boson arising 
in this model. 
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I. INTRODUCTION 


The CMS and ATLAS collaborations at the Large Hadron Collider (LHC) have suc¬ 
cessfully discovered a new resonance [1, 2] with a mass of 125 GeV [3], which has proper¬ 
ties consistent with the Higgs boson predicted by the Standard Model (SM). The signal 
strengths of this boson in its various final states are in good agreement with the SM 
expectations at la. The nominal variations in certain production and decay modes could 
be due to some physics beyond the standard model (BSM) or could simply be due to 
insufficient statistics. 

It is well known that the SM cannot be the final theory of nature. The successful 
explanation of the hierarchy problem requires some new physics (NP) near the TeV scale. 
In addition, the observation of small neutrino masses and their very particular mixing 
indicates the presence of physics beyond the standard model (BSM). There are a few well 
motivated theories, such as supersymmetric extensions of the SM or theories with extra 
spatial dimensions, that cure the aforementioned limitations of the SM. However, neither 
ATLAS nor CMS have yet conclusively discovered any particle that serves as proof for 
BSM physics. Now, with the discovery of the Higgs boson, effects of new physics can be 
searched for in its coupling measurements [4-37]. 

In this paper, we consider the simplest manifestation of a BSM extension through an 
extra singlet scalar. As a first step, we would like to see how the addition of just an 
additional neutral Higgs boson fares with the discovery prospects at the high-luminosity 
run at LHC (HL-LHC) with a final integrated luminosity of 3000 fb _1 . 

The presence of a heavy Higgs-like neutral scalar is innate in various models, such 
as, the minimal supersymmetric standard model (MSSM), two Higgs doublet models 
(2HDMs), models with extra spatial dimensions, etc. However, the simplest among these 
models is the SM augmented with a gauge singlet. This can originate very naturally 
from a U(1)b-l model with an extra 1/(1) local gauge symmetry [38], where B and L 
represents the baryon number and lepton number respectively. In particular, we focus 
on a TeV scale B — L model, that can further be embedded in a TeV scale Left-Right 
symmetric model [39-43]. The B — L symmetry group is a part of a Grand Unified Theory 
(GUT) as described by a S'O(IO) group [44], Besides, the B — L symmetry breaking scale 
is related to the masses of the heavy right-handed Majorana neutrinos, which participate 
in the celebrated seesaw mechanism [45-48] and generate the light neutrino masses. 

Another important theoretical motivation of this model is that the right handed neu- 
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trinos, that are an essential ingredient of this model participate in generating the baryon 
asymmetry of the universe via leptogenesis [49]. Hence, the B — L breaking scale is 
strongly linked to leptogenesis via sphaleron interactions that preserve B — L. ft is im¬ 
portant to note that in the U(1)b-l model, the symmetry breaking can take place at 
scales much lower than that of any GUT scale, e.g. the electroweak (EW) scale or TeV 
scale. Because the B + L symmetry is broken due to sphaleron interactions, baryogenesis 
or leptogenesis cannot occur above the B — L breaking scale. Hence, the B — L breaking 
around the TeV scale naturally implies TeV scale baryogenesis. 

The presence of heavy neutrinos, a TeV scale extra neutral gauge boson and an ad¬ 
ditional heavy neutral Higgs, makes the model phenomenologically rich, testable at the 
LHC as well as future e + e~ colliders [50-63]. The Majorana nature of the heavy neutrinos 
can be probed for example through same-sign dileptonic signatures at the LHC [64] . On 
the other hand, the extra gauge boson Z' in this model interacts with SM leptons and 
quarks. Non-observation of an excess in clilepton and di-jet signatures by ATLAS and 
CMS have placed stringent constraints on the Z' mass [65-70]. 

In this work, we examine in detail the discovery prospects of the second Higgs at the 
HL-LHC for a TeV scale U(1)b-l model. The vacuum expectation value (vev) of the 
gauge singlet Higgs breaks the U(1)b-l symmetry and generates the masses of the right 
handed neutrinos. We consider the B — L breaking scale to be of the order of a few 
TeVs, for which the right handed neutrino masses can naturally be in the TeV range. 
The physical second Higgs state mixes with the SM Higgs boson with a mixing angle 6, 
constrained by electroweak precision measurements from LEP [71-73], as well as from 
Higgs coupling measurements at LHC [74, 75]. The second Higgs is dominantly produced 
by gluon fusion with subsequent decay into heavy particles. The largest branching ratios 
are into W, Z and Higgs bosons. We discuss in detail the different channels through 
which the second Higgs state can be probed at the HL-LHC. 

Note that there are other possible B — L extensions of the SM, that have been studied 
in Refs. [76], Refs. [77] and Refs. [78]. In Refs. [76], the B — L gauge boson Z' acquires 
mass through the Stueckelberg mechanism [79, 80]. In this case, the B — L symmetry is 
unbroken, even after Z' acquires mass. Hence, the neutrinos in this model are necessarily 
of Dirac nature. To generate the mass of Z' via Stueckelberg mechanism, the presence of 
an axionic scalar is required. In addition to the U(1)b-l, an additional U(l)x symmetry 
is imposed, that brings down the scale of the Z l around TeV. As a second option [78], 
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the B — L symmetry is broken spontaneously by an SM gauge singlet Higgs field and as a 
consequence the Z' acquires mass. However, due to non-trivial B — L charge assignment 
of the gauge singlet Higgs, this scenario does not contain any Majorana mass term of the 
heavy right handed neutrinos. The light neutrinos in this model are again necessarily of 
Dirac nature. The collider signatures of these models are very different compared to the 
B — L extension where the light and heavy neutrinos are of Majorana nature. 

The paper is organised as follows: in section II, we review the basics of the U(1)b~l 
model. We discuss the constraints on the heavy neutrino sector and the limits on Z' 
in section III. Following this, in section IV, we outline the different constraints on the 
mixing angle between the SM-like Higgs and the second Higgs state. We study in detail 
the collider signatures of the heavy Higgs in section V. We briefly discuss non-standard 
production of the heavy Higgs in section VI. Decay of the heavy Higgs to a pair of heavy 
neutrinos is discussed in section VII. Eventually we offer conclusions in section VIII. 

II. BRIEF REVIEW OF THE U(\)b~l MODEL 

The U (1 )b-l model is one of the simplest extensions of the SM [44, 81-85]. In addition 
to the symmetry group of the SM, it has an additional U{ 1) gauge symmetry, that is 
identified as B — L symmetry. The full group structure of this model is therefore 


SU (3)c x SU(2) l x U (l)y x U{1) B -l, 


( 1 ) 


where U(1 )b-l represents the additional gauge symmetry. The Lagrangian of this model 
is as follows: 


( 2 ) 


C — C s + Cym + + £y, 


where C s , Cym, Cf and £y are the scalar, Yang-Mills, fermion and Yukawa terms respec¬ 
tively. The different terms in the Lagrangian are explained in detail in Refs. [51, 52, 57]. 
The Yang-Mills Lagrangian can be expressed as 



(3) 


where the first three terms represent the kinetic terms of the SU(3)c, SU{2 )l and U( l)y 
gauge groups respectively, a, b are the colour and SU (2) indices respectively. The fourth 
term is the kinetic term for the U(1)b-l gauge group and is represented by 



( 4 ) 
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where B' is the U(1) B -l field strength. 

In addition to the standard particle contents of SM, the fermion sector of this model 
has three right-handed neutrinos Nr, that are singlets under SM gauge group. This is 
required for anomaly cancellation and these right handed neutrinos generate Majorana 
masses of the light neutrinos through the seesaw mechanism, as discussed in section III. 
Analogous to the SM, the covariant derivative for this model is defined as 

D, = d, + ig s t a G + igT b W b + i 9l YB, + iejY B - L B^ (5) 

where g s , g, g\ and g' are the SU(3)c , SU(2)l , U(l)y and U(1)b-l couplings with t a , T b , Y 
and Y b _ l being their respective group generators. In the present study, we explicitly 
assume that there is no direct mixing between the two U( 1) fields B and B'. This 
corresponds to the minimal version of the B — L model. The fermion sector of the 
Lagrangian is expressed by 

£/ = {i (QL)iY D v(Q L )i + * { u R)il^ l D ll {uR)i + i (d^i^Dn(dR)i 

i= 1,2,3 

+2 + i (e^i^D^e^i + i ( NR) i 'y ,1 D fl (NR)i ), (6) 

where the electromagnetic charges on the fields are the same as the SM ones and the 
n — L charges are Y^-l = f anc ^ ^ b-l — — 1- 

In order to break the B — L gauge symmetry and to generate the mass of the additional 
Z' boson corresponding to this extra gauge symmetry, one needs to introduce an extra 
complex Higgs field y. The field y is a singlet under the SM gauge group with non-zero 
B — L charge Y£_ l = +2. The B — L symmetry is broken spontaneously by the vev of 
the Higgs field y. The SM Higgs field is neutral under the B — L gauge group, hence it 
has Y b _ l = 0. This particular choice preserves gauge invariance. 

The most general and renormalisable scalar Lagrangian of this model can be expressed 
as 


= (d^hYd^h + (D^xY(D,x) - V(x, H ), ( 7 ) 

where the scalar potential V (y, H) has the following form, 

V(x, H) = M 2 h H^H + ml |yf + A 1 (i7 t i7) 2 + A 2 |y| 4 + A 3 (^ t ^)|y| 2 . (8) 

To complete the discussion on the Lagrangian, we write down the Yukawa term, which 
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in addition to the SM terms has interactions involving the right-handed neutrinos Nr, 

c y = - - v% WuMiii - vtjiMMjH 

- N R )jH - y%Wrf<(N R ) jX + h.c., (9) 

where H = ia 2 H* and i,j runs from 1-3. The vev of the x field breaks the B — L 
symmetry and generates the Majorana masses of the right handed neutrinos Nr where 
M Nr = y M v'. On the other hand, the Dirac masses for the light neutrinos are governed 
by the Yukawa couplings y u s. 

Next, we turn our attention to spontaneous electroweak symmetry breaking (SSB) in 
this model. Further details are discussed in Ref. [57]. In order for the potential to be 
bounded from below, the couplings Ai ;2 ,3 should be related as 

4AiA 2 — A 2 > 0, 

Ai, 2 > 0. (10) 


In order to minimise the potential, one requires the above two conditions to hold. 

We denote the vevs of H and x by v an d v' respectively. On minimising the potential 
V(x,H) with respect to both vevs, one obtains [51, 52, 57], 

, 4A 2 M 2 - 2X 3 M 2 ;2 4AiM 2 - 2A 3 M 2 

v 2 = H _ v ' 2 = j H (11) 

A| — 4AxA 2 ’ Ag — 4 AxA 2 ' K } 

Since the B — L breaking scale is higher than the electroweak symmetry breaking 
(EWSB) scale, we have v' > v. 

X mixes with H due to the A 3 -term as shown in Eq. 8. The mass matrix between the 
two Higgs bosons in the basis (H, x) is given by 


M(H, X ) = 2 


\\v 2 X 3 vv'/2 

X 3 vv'/2 X 2 v' 2 


( 12 ) 


Next, we define the mass eigenstates as (Hi, H 2 ) which are related to the (H, x) basis by 



cos 6 — sin 6 
sin 6 cos 6 



(13) 


where the mixing angle 6 (—| < 6 < |) satisfies 


tan 26 


X 3 v'v 

(X 2 v' 2 — X\v 2 ) 


(14) 
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The masses of the physical Higgs bosons, Hi and H 2 are 


Mj Ii = Xiv 2 + X 2 v' 2 — \J (Ain 2 — X 2 v’ 2 ) 2 + A|n /2 n 2 , 
m h 2 = Xiv 2 + X 2 v' 2 + \J (Ain 2 - A 2 n' 2 ) 2 + X\v ,2 v 2 . (15) 

After imposing SSB in C s , one obtains the mass spectrum of the gauge bosons 


M 1 = 0, 

M w± = ^vg, 

M z = l ^\fg 2 + ~gl, 

Mz ' = 2 v'g'. (16) 


Note that, among the two Higgs masses, one chooses m 2 Hi < m 2 H , 2 , be., H\ is chosen 
to be the lightest state. In our subsequent discussion we will consider the case, where 
Hi is SM-like with a mass around 125 GeV. The other Higgs state, H 2 is heavy and we 
allow its mass to vary in the range 250-900 GeV for our phenomenological studies. The 
interactions of Hi and H 2 with the fermions and gauge bosons are expressed in terms of 

the mixing parameter 6 in the following manner 

eMf cos 0 


Hiff: 
HiW + W 


2M\y 
M\ye cos 6 


HiZZ : 


Syj 

M\ye cos 6 

9 ’ 

C W S W 
/2 / 


H 2 ff : 


H 2 W + W~ 


HoZZ 


HiZ'Z' : -8sin^'V, H 2 Z'Z' : 


eMf sin 6 
2M W 
Mwe sin 6 

Mwe sin 6 

9 ? 

C w S w 

-8 cos 0g' 2 v'. 


(17) 


The scalar self-interactions are given by 

HiHiHi : — 3-(4cos 3 9sin9 w M w Xi — 2sin 3 0eX 2 v'— 

e 

cos 2 9 sin 9eX 3 v' + 2 sin 6 W sin 2 9 cos 9M W A 3 ), 

H 2 H\ Hi : —(12 cos 2 9 sin 9 W sin 9MwXi + 6 sin 2 9 cos 9eX 2 v'+ 
e 

(1 — 3 sin 2 9) cos9eX 3 v' — 2(2 — 3 sin 2 9) sin^sin^M^/As). 


(18) 


In the above expressions, / denotes the SM fermions. We refer the readers to Refs. [51, 
52, 57] for a detailed description of the other interaction terms, arising from this model. 
Since, there are no extra coloured or electromagnetically charged states in this model that 
can alter the loop functions, we calculate the effective vertices ggH L 2 , lX'H i 2j Z^H X 2 and 
Z'^H\ f 2 following the standard loop functions relevant for the SM (see [86] and references 
therein). 
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III. CONSTRAINTS ON HEAVY NEUTRINOS AND Z' 


Before proceeding to discuss the phenomenological aspects of the heavy Higgs, we 
briefly discuss the constraints and limits on the various parameters arising from the 
heavy neutrinos Nr and Z'. 

A. Constraints on Heavy Neutrinos 

As we discussed in the previous section, the model consists of three right-handed 
neutrinos Nr, required for anomaly cancellation in the theory. The low-scale breaking 
of the B — L gauge symmetry implies the right handed neutrinos to be of the order of 
a few hundred GeVs to a few TeVs. In this work, we consider the heavy right handed 
neutrinos of TeV scale, that naturally emerge from low scale B — L breaking, without 
any unnatural tuning of the Yukawas. This scale will be accessible in the coming runs of 
the LHC and possibly also at future lepton and hadron collider experiments. The right 
handed neutrinos generate the masses for the light neutrinos via the seesaw mechanism 


M„ = -MjMpMc. 


(19) 


In the above, Mr> is the Dirac mass matrix of light neutrinos, whereas Mr is the Majorana 
mass matrix of the heavy neutrinos. By demanding Mr ~ TeV and M. u ~ eV, one is 


able to constrain the active-sterile neutrino mixing Mr/Mr ~ 10 6 . In addition to the 


constraints from light neutrino masses, the active-sterile mixing can also be constrained 
from other experimental searches, e.g. the neutrino-less double beta decay (Oz/2 f3), f3- 
decay, peak searches and kink searches [50, 63, 87]. The heavy Majorana neutrino, that 
mix with the active light neutrino v with a mixing angle 6 U , participate in the Oz/2/3-decay, 
where the amplitude is expressed as 



( 20 ) 


The non-observation of any positive signal in this lepton number violating process con¬ 
strains the active-sterile mixing to 6/ < 10 -5 , for a heavy neutrino mass Mr ~ 500 GeV 
[88]. A complete discussion on the different bounds on the active-sterile neutrino mixing 
can be found in [87, 89]. In addition, the collider signatures of the heavy neutrinos at 
LHC has been discussed in details in Ref. [54], A detailed discussion of the like sign 
dilepton signature from right handed neutrino decay has been studied in Ref. [54], 



B. Limits on Z' 


The B — L model has an additional gauge boson Z' of mass M' z = 2v'g[. Z’ interacts 
with the leptons, quarks, heavy neutrinos and light neutrinos with interaction strengths 
proportional to the B — L gauge coupling g’. The Z’ boson can in principle be detected by 
observing di-leptonic and di-jet signals at colliders. The presence of a sequential SM-like 
(SSM) Z l gauge boson has been severely constrained by direct searches at colliders, as 
well as by indirect searches. The ratio of Z’ mass to its coupling is constrained from 
indirect searches to be around [52, 77, 90, 91] 

717 yi 

— > 6.9 TeV. (21) 
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Several studies have been carried out by ATLAS and CMS in di-leptonic and di-jet 
channels to search for this elusive heavy gauge boson [65-70]. The Z' can decay to a 
boosted tt pair which provides sensitivity in semi-leptonic or fully hadronic top decays [69, 
92] . The cross-section times branching ratio ( crxB) has been constrained to be less than 1- 
2 pb [69] for a Z' with a width to mass ratio between T Z //M Z i = 1% and T z >/M Z i = 10%. 
Note that, here and in Table. I, we quote the most conservative limits of the cross-sections, 
where for other different masses the cross-sections are even more stringent. The recent 
combined analysis by CMS for the di-electron and di-muon mass spectra has further 
constrained the ratio ( R ) of cross-section times branching ratio of a narrow resonance [67]. 
In addition, the di-leptonic search by ATLAS has also constrained the sequential Z’ [66]. 
The other searches correspond to 

1. The search for a cli-jet resonance by CMS [68] that constrains a x B x A < 0.2 — 0.3 
pb (A being the acceptance for the kinematic requirements) and M z i <1.70 
TeV 

2. ATLAS search for t + t~ pair [70] 

3. CMS search for heavy resonance into bb pairs that bounds M Z ’ ssm [65]. 

In addition, we also show the limits applicable for a B — L model by comparing the 
limits from the 8 TeV run of the ATLAS di-lepton search [66], in Fig. 1. We consider few 
benchmark values for the free parameter g' and also for the mass of the heavy neutrino. 
Note that the production cross-sections of Z' in the B — L model have been computed 
at leading order (LO). The bounds on M Z i from this model are summarised in Table II. 
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(a) 


(b) 


FIG. 1. The comparison between the limits from ATLAS di-lepton search [66] with the B — L 
predictions with (a) Mjy = 500 GeV and g' = 0.05 (brown fine dotted), g' = 0.1 (magenta 
dashed) and g' = 0.2 (blue dotted) and (b) g' = 0.2 and Mjy = 500 GeV (blue dotted) and 
Mjv = 1 TeV (magenta dashed). 


In Fig. la, we find that by varying the coupling parameter r/, the bound on Mz> changes 
considerably by a few 100 GeV, whereas in Fig. lb, we find that by varying the masses 
of the heavy neutrinos, the bounds shift by (9(10) GeV 1 . Finally, we give an estimate of 
how the signal fares with respect to the SM backgrounds by studying the di-lepton final 
state using a basic set of trigger cuts on the transverse momentum ( p T ), pseudo-rapidity 
(rj) and isolation in the pseudo-rapidity-azimuthal angle plane (A R), i.e. pr,i > 10 GeV, 
|7^| < 2.5 and A Ru > 0.2. For the 14 TeV run, the benchmark Mz' = 3 TeV and 
g' = 0.2 yields the LO cross-section to be around 1.4 fb. Whereas, for the background 
the LO cross-section is around 1900 pb, several orders of magnitude larger than the signal 
cross-section. However, by imposing a simple invariant mass cut on the dilepton system, 
2900 GeV < Mu < 3100 GeV, one sees a dramatic reduction in the SM background, 
which amounts to ~ 0.01 fb. The signal, however reduces by a nominal amount to ~ 1.3 
fb. Hence, a massive Z' boson has a significant discovery potential during the 14 TeV 
LHC runs. 


1 A recent study [93] showed the importance of exclusion plots as functions of both masses and couplings. 
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Searches 

Constraints 

M z , (SSM) 

Boosted tt [69] 

a x B < 1 — 2pb 

- 

di-lepton-CMS [67] 

R < 7 x 10 -6 

2.90 TeV 

di-lepton-ATLAS [66] 

(j x B < 4 x 10~ 2 pb 

2.90 TeV 

di-jet-ATLAS [68] 

a x B x A < 0.2 — 0.3pb 

1.70 TeV 

r+T-- ATLAS [70] 

a x B < 0.1 pb 

1.90 TeV 

bb- CMS [65] 

- 

1.20- 1.68 TeV 


TABLE I. The recent bounds on Z' production from di-lepton, di-jet and other analyses. 


M n (TeV) 

g' 

M Z ,(B - L) (TeV) 

0.5 

0.2 

2.62 

1.0 

0.2 

2.65 

0.5 

0.1 

2.25 

0.5 

0.05 

1.83 


TABLE II. The bounds on Mz> derived from the ATLAS di-lepton search [66] relevant for a 
U(1) B -l model. 

IV. CONSTRAINTS ON HIGGS MIXING 

As discussed in the previous section, the light Higgs, H\ (or the SM-likc Higgs) and 
the heavy Higgs, H 2 mix with an angle 6. Hence, their couplings to the other particles 
in the model are scaled accordingly. Before discussing the phenomenological aspects of 
the searches, we impose bounds on the mixing parameter from the available experimental 
results. There are further theoretical bounds on this parameter which we discuss below. 

• Experimental bounds : Recent searches from CMS [74] and ATLAS [75] have 
already put bounds on a large class of BSM models. We work in the so-called n 
framework, where the coupling deviations of the SM-like Higgs are parametrized in 
terms of simple rescalings. The Higgs coupling to two fermions gthff and two weak 
bosons gH x vv are defined as [94], 

dHxff = Af-S'lr// and gn x vv = Kv-gmrvi ( 22 ) 

where htf and Ky are the coupling modifiers and are equal to unity in SM. We quote 
the 95% CL intervals on the various k parameters from CMS (Fig. 12 in Ref. [74]) 
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and ATLAS (Fig. 15 in Ref. [75])in Table III. Here, the experimental collaborations 
have assumed that the loop level couplings like Higg, H 177 and HiZ'y can be 
parametrized in terms of the tree level couplings and that no new loop particles are 
involved. They also assume that the invisible branching ratio of H\ is zero. These 
assumptions agree with the model under consideration. We consider the mass of 
heavy neutrinos to be in the TeV scale, for which the SM like Higgs decay to heavy 
neutrinos is small. In addition, since the light and heavy sterile mixing is small 
~ 10 -6 , hence, the SM like Higgs decaying to two light neutrinos is negligible. In 
this particular model, the couplings K t = Kb = Kw = Kz = k t = cos 9 and can have 
a maximum value of unity. Thus, all the major production cross-sections for Hi, 
e.g. ggF, VBF , VH and ttH scale as cos 2 6. 


KW 

KZ 

Kb 

k t 

CMS 

[0.66,1.24] 

[0.69,1.37] 

[0.51,1.22] [0.07,1.46] 

[0.47,1.25] 

ATLAS 

[0.63,1.19] [ 

-1.20, -0.67] U[0.67,1.26] [0.59,1.39] [-1.29,1.31] [ 

-1.46, -0.61] U[0.62.1.47] 


TABLE III. The 95% CL ranges on various signal strength modifiers, k, as reported by CMS [74] 
and ATLAS [75]. 


Using the ranges in Table III, we obtain the scale factor of the heavy Higgs, H 2 as 
sin 2 9 < 0.31(0.33) for CMS (ATLAS) at 95% CL. 

It is however important to note that the bounds on sin 6 from the coupling mea¬ 
surements of the SM-like Higgs are possibly the most desired and robust ones. 
These bounds are independent of the mass of the heavy higgs and will probably 
get more stringent with more integrated luminosity. As an example, in Ref. [95], 
the H —y WW* measurement is shown to constrain sin 6 ~ 0.36 from the projected 
study of LHC at 14 TeV with J Cdt =300 fb _1 . The same study also projects a 
smaller sin 9 ~ 0.25 at the ILC, running at 250 GeV with an integrated luminosity 
of 250 fb _1 . The ILC runs with greater centre-of-mass energies and higher inte¬ 
grated luminosities are expected to constrain sin 6 to even smaller values. In this 
analysis, we assume sin 6 = 0.2 which is in sync with the projected study at LHC 
14 with 300 fb -1 . 
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• Theoretical bounds : 


— Constraints from Mw • One of the strongest constraints on the the mixing 
angle, sin0, comes from the one-loop correction to the IV-boson mass, Mw , 
which is required to agree within 2cr of its experimental value, i.e., Mw = 
80.385 ±0.015 GeV [73, 96, 97]. This has recently been studied in the context 
of this model in Refs. [71, 72, 98]. It has also been shown in Ref. [71], that in 
the high mass region, the constraints from the one-loop correction to Mw are 
stronger than the ones obtained from S,T and U parameters [99-102], The 
upper bound on sin# decreases from ~ 0.35 to ~ 0.20 as Mh 2 increases from 
250 GeV to 900 GeV [71]. In our analysis, we have considered a conservative 
value of sin# = 0.20, throughout, in order to satisfy all the constraints. 

— Constraints from perturbative unitarity : Demanding perturbative uni- 
tarity [103], by studying all the 2 —)• 2 scattering amplitudes and demanding 
that the partial wave amplitudes a^s follow 

|R*M < (23) 

where the subscript £ denotes the orbital angular momentum, results in an up¬ 
per bound on the Higgs boson mass. The bounds from perturbative unitarity 
for a model with a scalar extension has been derived in Ref. [104], Perturbative 
unitarity also poses strong constraints on the ratio tan/3 = v/v'. 

— Perturbativity of the couplings : All the couplings in the potential are 
required to conform within perturbative limits, i.e., Ai t 2,3 < 47T. These bounds 
are weaker than the ones obtained from perturbative unitarity at the EW scale. 
Besides, constraints from vacuum stability and the renormalisation group evo¬ 
lution of Ai j2 ,3 are also studied in Refs. [62, 71, 105]. 

V. COLLIDER SEARCHES FOR THE HEAVY HIGGS 

The heavy Higgs H 2 in the B — L model mixes with the SM-like Higgs, H 1 , with 
mixing angle 6 , as has been discussed in the previous sections. H 2 can be produced at 
the LHC through multiple production processes, e.g. gluon fusion ( ggF ), weak boson 
fusion ( VBF ), associated WH 2 /ZH 2 productions and the associated ttH 2 production 
mode. Once produced, H 2 promptly decays into different final states, with WW, H 1 H 1 
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and ZZ being the dominant decay modes. In this section, we study in detail the collider 
signatures of H 2 produced through its dominant production mode, ggF , after including 
the constraints on the mixing angle, 9 , as discussed above. In order to study the col¬ 
lider signatures of H 2) we implement the model using FeynRules [106]. The generated 
Universal FeynRules Output (UFO) [107] hies are then fed to the Monte-Carlo (MC) 
event generator MadGraph [108] for generation of event samples. The parton-showering 
and hadronisation is carried out in the Pythia 6 [109] framework. For jet formation, we 
use the anti-fc^ algorithm with a jet parameter of R = 0.4 [110]. 

In Fig. 2a we show the branching ratios of H 2 to various final states as function of 
its mass, varying 1W # 2 between 250 GeV and 1 TeV. As is clear from Fig. 2a, the three 
most dominant decay modes of H 2 are WW, HiHi and ZZ. In Fig. 2b, we show the 
Next-to-Next-to Leading Order (NNLO) cross-sections of the three different final states 
mentioned above. Note that, in addition to the aforementioned processes, we also show 
the cross-section for the process pp —>■ H 2 —$■ WW —> 2£$ T in Fig. 2b. However, we 
do not consider the phenomenology for the latter process because of a somewhat less 
amount of handle on its kinematics due to fewer visible particles in the final state. The 
cross-section of H 2 decaying to lv2j is the highest, whereas for the 4£ channel, the cross- 
section is the smallest. We analyze these two channels in considerable details and study 
the discovery prospects of H 2 at the HL-LHC. We also briefly mention the 2t2j final state 
as a potential channel for discovering H 2 . 

Recently, search strategies for H 2 —>■ H\H\ have been discussed in Refs. [98, 111-114]. 
CMS [115] and ATLAS [116, 117] have studied the di-higgs production in the bbbb and 
6677 final states mostly in the context of models with extra spatial dimensions. The upper 
limits on a x B for the resonant and non-resonant production of di-higgs in context of 
such BSM models are found in Refs. [115-117]. A naive leading order estimate of the 
pp -7- H\H\ cross-sections [118] with v' = 3.75 TeV and sin$ = 0.2 reveals that for lower 
values of Mh 2 , he., up to ~ 500 GeV, the cross-section is substantially enhanced with 
respect to the SM cross-section. However, with higher values of Mh 2 , H 2 decouples and 
the cross-section tends to the SM value, see Fig. 3. Hence, this channel can complement 
the gauge-boson final states in searches for H 2 . 

In the present work, we focus on the WW and ZZ decay modes and try to devise 
some search strategies in the context of the 14 TeV run at the HL-LHC. Although, the 
branching ratios oiW/Z decaying to di-jet final states are large, still, the leptonic and 
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(a) (b) 

FIG. 2. Left panel : Branching ratios of H 2 (as a function of Mh 2 ) to W + W~ (solid magenta), 
ZZ (dotted red), H\H\ (dot-dashed black), tt (dotted black), bb (solid black), gg (solid red), 
rf (dotted blue), Z 7 (solid green), cc (fine-dotted red), 77 (solid blue), pp (fine-dotted green) 
and ss (fine-dotted black). Right panel : NNLO Cross section (fb) times Branching ratio as 
functions of Mh 2 (pp H 2 —> W + W~ —> Ivjj [solid red], pp -> H 2 —> W + W~ —> £ + v£~u 
[solid blue], pp —> H 2 —> ZZ —> £ + £~jj [solid green] and pp —> H 2 -> ZZ —> M [dotted black]). 
sin@ = 0.2 for all the cases. 



FIG. 3. LO cross-section for the channel pp — > H\H\ for (a) the U(1)b~l model (solid red 
curve) and (b) the SM with Mf l = 125 GeV (green curve). 
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semi-leptonic decay modes offer the cleanest possible signatures because of significantly 
less backgrounds. Therefore, in our subsequent discussions, we concentrate only on those 
channels, that have leptons in the final state, i.e. 

• pp —» H 2 —» ZZ —> 4£, 

• pp —>■ H 2 ZZ 2j2P and 

• pp —y H 2 —> WW —> lu2j 

For the search strategies, we adopt two different reconstruction methods which we 
discuss below. 

• Cut-based analysis (CBA) : In this method, we employ rectangular cuts on 
various kinematic variables in order to optimise the significance 

n = Afs / \/A5 + A Tq ■ 

• Multivariate analysis (MVA) : We employ multivariate techniques for better 
signal-to-background discrimination, resulting in better signal significance, n. For 
the present study, we use the Boosted Decision Tree (BDT) algorithm from the 
TMVA [119] framework. In order to perform an MVA, we select the set of kinematic 
variables that give the maximum discrimination between signal and background. 
Both the signal and backgrounds are trained by this algorithm and another set 
of event samples are used to test the BDT output. For any MVA, we must al¬ 
ways be alert not to over-train signal and background. The universally accepted 
Kolmogorov-Smirnov (KS) test can reveal if our choice of parameters needs to be 
changed. The test sample is not over-trained if the KS probability lies in the 
range (0.1, 0.9). For most cases, a critical value of the KS probability greater than 
0.01 [120] implies that the samples are not over-trained. For the subsequent studies 
we ensure that over-training is not an issue over the entire parameter range. In 
order to estimate the LHC’s potential in excluding H 2 , we use the result of the 
MVA as input to a binned log-likelihood hypothesis test [121]. 

In the following subsections we discuss the discovery prospects of the heavy higgs. 
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A. pp -> H-2 ->• ZZ -> 4f 

In this scenario, H 2 is produced on-shell and decays to two Z bosons. The two Z 
bosons subsequently decay to four leptons. The main background for this process is the 
ZZ production mode that will generate the same final state. To analyse this channel we 
employ the following trigger cuts: 

• Trigger Cuts (TC): 

To identify the leptons, we apply the following minimal cuts. 

1. Transverse Momentum: px(l) > 10 GeV 

2. Pseudo-rapidity: \rj(l)\ < 2.5 

3. Radial Distance: A R(li,lj) > 0.2 

We show the normalised distributions of various kinematic variables in Figs. 4 and 5 
for M H2 = 250, 500 and 900 GeV with respect to the background. It is evident 
from Fig. 5a, that while for the signal, the invariant mass of four leptons peaks at 
the heavy Higgs mass, Mh 21 this is not true for the background and hence serves 
as one of the better discriminating variables. The leptons originating from the Z 
decays also have higher transverse momentum compared to the background, which 
peaks at lower px values. The two Z bosons also have higher px with respect to 
the background, as shown in Figs. 5b and 5c. For higher M H . 2 , both the Z s have 
Pt > 100 GeV. For our cut-based analysis (CBA), we use these following selection 
cuts to separate signal from background to a good degree. 

• Selection Cuts (SC) 

We use the following selection cuts: 

1. Invariant mass of the four lepton system: to lie in the range, Mh 2 ± 10 

GeV 

2. Transverse momentum of leading lepton: px ti > 90 GeV 

3. Transverse momentum of sub-leading lepton: px e > 70 GeV 

4. Transverse momentum of the other two leptons: px tz > 50 and pr l/L > 20 GeV 

5. Invariant mass of the reconstructed Z bosons: Mz ,, Mz 2 E Mz ± 10 GeV 

6. Transverse momentum of the two reconstructed Z bosons: Pt{Z\), Pt(Z 2 ) > 
100 GeV 
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FIG. 4. pp —> ZZ —> M channel normalised distributions for Mh 2 = 250 GeV (green dotted 
line), 500 GeV (blue dashed line), 900 GeV (red solid line) and background (purple solid): 
(a)-(d) pt distributions of the four leptons, pt sorted. 


For the MVA analysis, we choose a set of 18 kinematic variables with the maximum 
discriminating power, which are M 4 ^, pr, , A Mz k , px(Zk), rj(Zk) and pr(M). Here 
i, j = 1,4, fc = 1,2 and the leptons and two Z s are pr sorted. pr(4£) is the vector sum 
of p T of the four leptons. 

We tabulate the signal and background cross-sections after the trigger cuts (TC) and 
the selection cuts (SC) in Table. IV for Mh 2 in the range 300 — 700 GeV. We also list the 
number of signal (Ms) and background (Mb) events computed for an integrated luminosity 
of 3000 fW 1 for the cut-based and BDT analyses after imposing the respective cuts. Fig. 6 
shows the normalised distributions for the signal and the background against the BDT 
response for two benchmark masses, M# 2 = 250 GeV and 500 GeV. We find that with 
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FIG. 5. pp —> ZZ —> M channel normalised distributions for Mh 2 = 250 GeV (green dotted 
line), 500 GeV (blue dashed line), 900 GeV (red solid line) and background (purple solid): 
(a) invariant mass of the four leptons, (b)-(c) transverse momenta of the two reconstructed Z 
bosons and (d) vector sum px of the four lepton system. 


an increase in mass, the overlap between the signal and the background decreases. As a 
result, using a BDT, S/B improves significantly. 

The maximum significance is obtained for relatively small masses of II), where the 
cross-sections are sufficiently large. Note that, even with the cut-based analysis, the 
discovery prospect of the heavy Higgs is rich for 3000 fb -1 integrated luminosity. As 
an example, with the cut-based analysis, a Higgs with mass Mh 2 = 400 GeV can be 
discovered with 8.6cr significance {jicba) at the HL-LHC with an integrated luminosity 
of 3000 fb -1 . As expected, the BDT analysis is seen to improve the significance ( ubdt ) 
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FIG. 6. Normalised signal and background distributions against BDT response for (a) Mjj 2 = 
250 GeV and (b) Mh 2 = 500 for the channel pp —> H 2 —> ZZ —> A£. 

by a considerable amount. In the entire mass spectrum, the maximum difference in Ucba 
and ubdt occurs for M H2 = 300 GeV, where the signal and background distributions 
mostly overlap, making it very difficult to impose rectangular cuts. 

To estimate the necessary integrated luminosity to exclude the existence of the heavy 
Higgs boson, we use the BDT output, shown in Fig. 6 , weighted with the according cross 
section as input for a CLs likelihood ratio [121], see Fig. 7. Conservatively, we assume a 
flat systematic uncertainty of 10% for each bin. While an H 2 with Mh 2 = 250 GeV can 
be excluded at 95% CL with 100 fb _1 in this channel, excluding M H2 = 700 GeV requires 
3000 fb -1 . 

B. pp H 2 ZZ —> ‘It + 2 j 

The channel H 2 —)• ZZ —>■ 2/2 j has been studied in [122] in the context of heavy SM 
Higgs boson searches. While the signal benefits from a larger branching ratio of the Z 
boson to jets compared to leptons, the only major background in this channel remains 
continuum ZZ production. Here, for convenience, we briefly describe the selection cuts 
discussed in full detail in Ref. [122], 

• Leptonic Z reconstruction : We demand two isolated muons with px > 15 GeV 
and |? 7 | < 2.5. We further demand an invariant mass window of 10 GeV around 
M z . 
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FIG. 7. Confidence level contours for Mh 2 = (a) 250 GeV, (b) 500 GeV and (c) 700 GeV. 
We show results for integrated luminosities (f Cdt ) from 50 to 3000 fb _1 . We assume a flat 
systematic uncertainty on the backgrounds of 10%. 

• Hadronic Z reconstruction : We reconstruct the hadronic Z following the al¬ 
gorithm given in Ref. [123]. Here also we require an invariant mass window of 10 
GeV around M z . 

• Heavy Higgs reconstruction : If the previous two steps are successfully satisfied, 
then the invariant mass peaks as Mjj 2 = ( pz lep +Pz had ) 2 , where pz lep and pz had are 
respectively the four-momenta of the reconstructed Z bosons in the lcptonic and 
hadronic channels. The Higgs mass windows used for the four benchmark masses 
are (300 ± 30,350 ± 50,400 ± 50,500 ± 70,600 ± 100) GeV. These are found to 
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Mh 2 

(GeV) 

a TC 

(fb) 

a sc 

(fb) 

Mg BA 

M£ BA 

ncBA 

M BDT 

M BDT 

n BDT 

300 

0.126 

0.010 

30 

105 

2.62 

227 

555 

8.12 

350 

0.132 

0.042 

125 

162 

7.37 

262 

419 

10.03 

400 

0.113 

0.047 

142 

131 

8.60 

246 

361 

9.99 

450 

0.078 

0.034 

101 

101 

7.14 

168 

243 

8.29 

500 

0.051 

0.021 

63 

81 

5.26 

93 

132 

6.19 

550 

0.034 

0.013 

40 

48 

4.23 

54 

70 

4.82 

600 

0.022 

0.008 

24 

45 

2.87 

42 

112 

3.42 

650 

0.015 

0.005 

14 

32 

2.12 

23 

60 

2.54 

700 

0.010 

0.003 

9 

24 

1.57 

16 

87 

1.58 

SM 

28.626 

- 

- 

- 

- 





TABLE IV. NNLO cross sections after trigger cuts (cr ic ) and selection cuts (cr 5 ^). A fs and 
Mb represent the number of signal and background events, respectively, while the superscript 
and subscripts CBA and BDT represent the cut-based and BDT analysis, n is the significance. 
The number of events have been computed for an integrated luminosity 3000 fb _1 . All the 
cross-sections include the higher order corrections to the NNLO level. 

optimise the results. Also because the Higgs width increases significantly with 
Mh 2 , we widen the windows for reconstruction purposes. 

• ZZ separation : To further improve S/B, we require the leptonic and hadronic Z 
bosons to be have a maximum isolation of A Rzz < 3.2, where A R = A r\ 2 + A <j) 2 
with A rj and A0 being the separation between two objects in the pseudo-rapidity 
and azimuthal angle planes respectively. For Z + jets, A R between the recon¬ 
structed Zi ep and the fak e-Z from the QCD jets often becomes large in order to 
account for the large Higgs invariant mass. 

• Pruning and trimming : The pruning [124, 125] and trimming [126] algorithms 
are used to further reduce the QCD backgrounds because this technique helps in 
discriminating colour singlet resonances from QCD jets [127]. The details of this 
procedure are elucidated in Ref. [122], 

After assuming sind = 0.2 and applying the reconstruction outlined in [122] we find the 
results shown in Tab. V. 
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Mh 2 

(GeV) 

ggF+VBF 

a sc 

(fb) 

bkq 

a sc 

(fb) 

S/B 



S/VS + B wo 

S/y/S + -B3000 

300 

0.048 

2.10 

0.023 

0.331 

1.811 

400 

0.290 

19.21 

0.015 

0.657 

3.602 

500 

0.223 

18.01 

0.012 

0.522 

2.858 

600 

0.121 

11.83 

0.010 

0.351 

1.920 


TABLE V. crg 9 ^^ BF is the production cross-section of H 2 from the ggF and VBF channels 
combined after employing the selection cuts discussed in Ref. [122], is the background 

cross-section for the same set of selection cuts. The table also shows the discovery potential of 
H 2 in this channel with the help of S/B and S/y/S + B. The subscripts 100 and 3000 imply 
the significance computed at the respective integrated luminosities in fb -1 . 

Hence, the sensitivity in the H 2 —>• 2 £2j channel alone is fairly small for the U(1)b-l 
model, based on the reconstruction of boosted Z bosons. However, this channel can be 
combined with the other channels in a global fit. 


C. pp —> H 2 —¥ WW -> i + $ T + > 2 j 

In this scenario, H 2 decays to W + and W~, followed by the subsequent decay of one 
of the Ws to a lepton and missing energy and the other one to a pair of jets. 

Recently, the ATLAS collaboration has searched for WW/WZ resonances decaying to 
a lepton, neutrino and jets [128], that mimic our signal. With a pr(W) > 400 GeV cut, 
used in Ref. [128], our signal cross-sections will be extremely small. In Ref. [128], the 
authors have considered two benchmark models, viz., (i) a spin 2 Kaluza-Klein Graviton 
for the WW resonance and (ii) a spin-1 SSM W' decaying to a WZ pair. From Fig. 2b, we 
find that the production cross-section of pp —>■ H 2 —>■ W + W~ is approximately a < 0.098 
pb, that is way below the exclusion limit, as given in Fig. 2 of Ref. [128]. Hence, we adopt 
different sets of cuts which are suitable for our analysis. 

For a heavy Higgs boson, the intermediate Ws are expected to have large px- Because 
the W bosons are highly boosted, the leptons and jets are expected to have a large 
isolation of A R(l, j/). We show the normalised distributions of various kinematic variables 
in Figs. 8, 9, 10 and 11, for the masses M H2 = 250, 500, 900 GeV and for the SM 
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(d) 


FIG. 8. pp —> WW —>• 12j$T channel normalised distributions for Mh 2 = 250 GeV (green 
dotted line), 500 GeV (blue dashed line), 900 GeV (red solid line) and background (purple 


solid): (a) Invariant mass of l2j$ T i (b) px of lepton, (c)-(d) px of the two tagged jets. 


background. 

To identify the leptons and jets, we apply the following minimal cuts: 

• Trigger Cuts (TC) 

1. Transverse Momentum of the jets : px(ji) > 30 GeV, where i is the jet-index 

2. Transverse Momentum of the lepton : px(l) > 20 GeV 

3. Pseudo-rapidity of the lepton : \rj(l)\ < 2.5 

4. Pseudo-rapidity of the jets : |? 7 (jj)| < 5.0 

5. Radial Distance between jets i and j : ARij^jj) > 0.4 
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FIG. 9. pp —>• WW —>• 12j$T channel normalised distributions for Mh 2 = 250 GeV (green 
dotted line), 500 GeV (blue dashed line), 900 GeV (red solid line) and background (purple 
solid): (a) pt of the system, (b) px of the di-jet system which reconstructs the W mass, (c) 
Invariant mass of di-jet system which reconstructs the W mass, (d) scalar sum px of lepton, 
two-tagged jets and $x- 


6. Radial Distance between lepton and the i th jet: AR(l,ji) > 2.0 


In the above, j\ and j ’2 denote two leading jets, sorted according to their transverse 
momentum. For the Higgs masses, Mh 2 of our interest, both the intermediate gauge 
bosons are on-shell and this allows us to fully reconstruct them. For the hadronically 
decaying W, we use the jet four-momentum, while for the leptonically decaying W, we fix 
the longitudinal component of the neutrino momentum, p z (u) by imposing the constraint 
Mw = {Pi+Pvf- 
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FIG. 10. pp —» WW —> £2j$T channel normalised distributions for Mh 2 = 250 GeV (green 
dotted line), 500 GeV (blue dashed line), 900 GeV (red solid line) and background (purple 
solid): (a) Abetween the two tagged jets, (b) $ T , (c) pseudorapidity of the hardest px 
jet, (d) pseudorapidity of the second hardest px jet and (e) pseudorapidity of the lepton. 


The major part of the background originates from non-resonant W + W~ production, 
with one W decaying hadronically and the other decaying leptonically. From the differ¬ 
ent px distributions in Figs. 8 and 9, it is evident that for low masses, the signal and 
background has large overlap, making the discrimination a difficult task. In addition, for 
the invariant mass, the signal and background show a large overlap for Mh 2 ~ 250 

GeV (see Fig. 8a), whereas for larger masses, the overlap decreases, px of the lepton and 
the leading-jet are large for 500 GeV < M H2 < 900 GeV. We also show the transverse 
momentum of the two reconstructed W s in Figs. 9a and 9b. For the signal, they peak at 
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FIG. 11. pp —> WW —> £2j$T channel normalised distributions for M= 250 GeV (green 
dotted line), 500 GeV (blue dashed line), 900 GeV (red solid line) and background (purple 
solid): pseudorapidity of the lepton. 

Pt(W) > 100 GeV, while for the background, most of the intermediate Ws have lower 
Pt- In addition, the signal has a large missing transverse energy ($ T ), as evident from 
Fig. 10b. In Fig. 9c, we show the invariant mass distribution of the two hardest jets and 
in Fig. 10a we show the distribution of their azimuthal angle separation. 

Note that the partonic signal cross-section, <7 gig varies from few tens of fb to 0(0.1) 
fb, for M H2 varying between 300 GeV and 900 GeV. However, the background for this 
process is extremely large <7bkg ~ 3380 pb. Hence, to extract the signal from background 
in a statistically viable fashion, we categorise the signals into four separate regions and 
implement stringent cuts both at the generation level as well as at the detector level. 


Mh- 2 

(GeV) 

Pril/ji/h) 

(GeV) 

J2)min 

J2)max 

(GeV) 

350 

30 

0.4 

1.4 

50 

500 

40 

0.2 

1.0 

60 

700 

50 

0.2 

0.8 

70 

900 

70 

0.2 

0.6 

90 


TABLE VI. Basic trigger cuts used to separate the signal from background in the pp —>■ #2 —>■ 
WW —> £ + > 2 j channel. 

• The low mass region: M# 2 ~ 350 GeV case: For this case the signal and background 
has a large overlapping region. The transverse momentum variable for the lepton 
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Mh 2 

(GeV) 

PT,l/jl,2 

(GeV) 

PT,W lt2 

(GeV) 

a pmax 
AR h,h 

(GeV) 

S T 

(GeV) 

1 M ljj$ T ~ 

(GeV) 

\Mjj - M w | 

(GeV) 

350 

35 

100 

1.35 

55 

225 

50 

20 

500 

45 

100 

0.9 

70 

250 

50 

20 

700 

55 

100 

0.75 

75 

250 

50 

20 

900 

75 

100 

0.58 

95 

600 

50 

20 


TABLE VII. Selection cuts to separate out signal from the background in the pp —> H 2 —> 
WW —» £ + $t+ > 2 j channel. 
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* IKolmogorov-Smirnov test: signal (background) probability = 0.045 (0.711) 
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1: signal (background) probability = 0.041 (0.257) 
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(a) (b) 

FIG. 12. Normalised signal and background distributions against BDT response for (a) Mh 2 = 
350 GeV and (b) Mh 2 = 500 for the channel pp —> H 2 —> WW —> h'jj. 


and jets peak around 30 GeV for the background. Hence, to reduce the background, 
we implement selection cuts, such as, Pt > 35 GeV, missing transverse energy 
$ T > 55 GeV and the scalar sum of transverse momentum of the lujj system, 
St > 225 GeV. I 11 addition, the selection cuts on the invariant mass of the lujj 
system, |— Mh 2 \ = 50 GeV reduces the background significantly. 


• Intermediate mass range: M# 2 ~ 500 GeV: For the intermediate mass range, such 
as, Mh 2 = 500 GeV, the signal and background has a relatively smaller overlap, 
as is evident from Figs. 8b, 8c and 8d. In addition, the reconstructed Ws have 
large transverse momenta, as shown in Figs 9a and 9b. The invariant mass of the 
lujj system peaks around 500 GeV, while for background the invariant mass peaks 
at much lower values. The larger cuts on the transverse momenta of the Ws and 
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Mh 2 

(GeV) 

(fb) 

•^bkg 

(Pb) 

<4° 

(fb) 

CT SC 

(fb) 

^Z A « 

U CBA 

03000 

350 

16.85 24.22 5.49 1666.29 

549 

166629 

16474 

4998880 

500 

9.06 

6.69 

3.44 

360.26 

344 

36026 

10305 

1080769 

700 

2.39 

2.56 

0.84 

77.63 

84 

7763 

2509 

232887 

900 

0.57 

0.65 

0.06 

6.19 

6 

619 

193 

18564 


TABLE VIII. The NNLO signal (a s ) and background (obkg) cross-sections and the number of 
events for different Higgs masses. The ■^’sZZoj b 100 ' ^s-mao / B- mm are the signal and background 
events with 100 and 3000fb -1 integrated luminosities for the cut based analysis. 

the cuts on the invariant mass efficiently reduces the background cross-section from 
~ 3380 pb to ~ 6.69 pb. 

• Large masses: M# 2 ~ 700 GeV: For large masses, such as, Mh 2 = 700 GeV, the 
signal cross-section after the selection cuts is relatively small a ~ 0.84 fb. However, 
the background is well separated. Hence, the stringent cuts on the kinematic vari¬ 
ables improves the sensitivity. We use larger cuts on the transverse momenta of the 
lepton, jets and the reconstructed Ws, as shown in Tables VI and VII. The back¬ 
ground cross-section after the different trigger and selection cuts reduce to ~ 2.56 
pb. 

• Very large masses: M# 2 ~ 900 GeV: For Mh 2 as large as 900 GeV, the signal cross- 
section after imposing the selection cuts becomes extremely low, cr ~ 0.06 fb. Here, 
we thus use higher trigger and selection cuts as is given in Table. VI and Table. VII, 
in order to reduce the background. 

In order to perform the multivariate analysis, we choose a set of 27 kinematic vari¬ 
ables with excellent discriminating power, which are Pt(£), r/(£), Pr(ji), v(ji), $Ti 

</>(# T ), pr{e,$ T ), Pr{j 1 J 2 ), |A0(Wi,W 2 )|, |A0(£,j 1 )|, A rj(£,j 2 ), A r)(t,ji), |A0(j 1; j 2 )\, 
Ar/(ji, j 2 ), \A(j)(ji,$ T )\, S T , M hh M jlhh Aand AR(j 1 j 2 ). In the above, i = 1,2 
and the jets and the reconstructed W s are pr sorted. 

The results for the cut-based and the multivariate analyses are shown in Tables VIII 
and IX. As expected, we find the MVA to be performing better with respect to the cut- 
based analysis. Fig. 12 shows the normalised distributions of the signal and background 
as a function of the BDT response for two benchmark masses, i.e. M^ 2 = 350 GeV 
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Mh 2 


KZ 


(GeV) 


350 

591 

116793 

17731 

3503792 

500 

309 

18761 

9270 

562838 

700 

68 

3718 

2055 

111544 

900 

12 

1358 

366 

40735 


TABLE IX. The signal and background events after BDT cut for different Higgs masses. The 

are the signal and background events with 100 and 3000 fb 1 integrated 

‘-’100/i>100 ’ ‘->3000/£>3000 ° ° 0 

luminosities. 


and 500 GeV. We see that for both cases, there is a large region of overlap between the 
signal and background. This implies that even with a BDT, the separation of signal 
and background is challenging, resulting in a very small S/B < 1/100. Assuming zero 
systematic uncertainties we quote the statistical significance in Table X. The discovery 
potential of H 2 in this channel is however bleak. 


M H2 GeV 

C [fb’ 1 ] 

riCBA 

n BDT 

350 

100 

1.34 

1.73 

3000 

7.36 

9.45 

500 

100 

1.80 

2.24 

3000 

9.86 

12.26 

700 

100 

0.94 

1.11 

3000 

5.17 

6.10 

900 

100 

0.26 

0.33 

3000 

1.41 

1.81 


TABLE X. The significance for cut-based and multivariate analysis for integrated luminosity 
100 fb’ 1 and 3000 fb’ 1 . 


VI. NON-STANDARD HEAVY HIGGS PRODUCTION 

In addition to the standard production processes, i.e. gluon fusion, VBF , VH 2 (V = 
W, Z ), ttH 2 , the heavy Higgs H 2 can also be produced in association with a Z' [52, 57] 
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FIG. 13. Left panel: The LO cross-section for the associated production pp Z'H -2 for mixing 
0 = 0 and different values of Mz> and g' (solid red, dashed green, dotted blue and dot-dashed 
magenta for [Mz> = 1 TeV, g' = 0.145], [Mz> = 2 TeV, g' = 0.290], [Mz> = 3 TeV, g' = 0.435] 
and [ Mz■> = 2 TeV, g’ = 0.2] respectively) such that ^r- >6.9 TeV (Eq. 21). Right panel: 
Comparison between the associated production pp —> Z'H 2 between sin0 = 0 and sin0 = 0.2 
(dashed green and dotted blue are for [Mz> = 2 TeV, g' = 0.290, sin0 = 0] and [Mz> = 2 TeV, 
g l = 0.290, sin 6 = 0.2] respectively). 


J . In the decoupling limit, the mixing 9 ~ 0 and hence the gluon fusion contribution 
would be negligible. However, the vertex factor Z l Z l H 2 is proportional to cos 6 (see 
Eq. 17). Therefore, even in the decoupling limit, the process pp —> Z’H 2 will give a 
non-zero contribution. We show the production cross-section of this process in Fig. 13a 
for the decoupling limit. Note that, for a non-zero value of 9, both the s-channel ( 
pp — > Z’* —f H 2 Z' ) 3 and t-channel diagrams mediated by quarks will contribute (see 
Fig. 13b). However, the t-channel contribution will be small. 


VII. DECAY TO HEAVY NEUTRINOS 

Here we briefly discuss the decay of H 2 to two heavy neutrino states, i.e., pp —> H 2 —■> 
NrNr. The unique feature of the B — L model is that both Z' and H 2 can decay to 

2 For a complete list of production processes, see [51, 57, 58]. 

3 Note that, Z'H 2 r y contribution is negligible because the vertex is loop suppressed and also because in 
the minimal B — L scenario, the Z'W + W~ coupling is absent. Hence, we neglect this contribution in 
our study. 
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FIG. 14. Left panel: The NNLO production cross-section of pp — > H 2 —> NrNr. Right panel: 
The branching ratios of Ho —> NrNr (solid blue : Mr = 200 GeV and dotted red : Mr = 100 
GeV). 


a pair of heavy neutrinos. These produced heavy neutrinos can further decay to IW, 
vZ and H\v states, producing same as well as opposite sign leptonic signatures. The 
same sign leptonic signature confirms the Majorana nature of the heavy neutrinos, as 
well as, the presence of B — L breaking 4 . In Ref. [54] the authors have studied the 
same sign dilepton+4-jets signature from pp —> Z' —» NrNr channel. I 11 Ref. [55], 
multilepton final states have been studied. In Fig. 14a we show the NNLO cross-section 
of pp —* H 2 —>■ NrNr, which is only few fb for Mn = 100 GeV. In addition, we also 
show the branching ratio of H 2 — > NrNr in Fig. 14b. Note that, this branching ratio is 
small. Hence, even for lower Nr masses, e.g. 100 or 200 GeV, the previous analysis of 
pp —)■ H 2 —>■ W + W~ and pp —>■ H 2 —)■ ZZ will remain practically unchanged. A detailed 
study for pp —» h 2 —$■ NrNr has been presented in Refs. [57, 58]. This channel also offers 
other final states, e.g. + 4 j, / ± Z =F i ± + 2 j + $ T , 4 1 + $ T , 4b + $ T and l + $ T + bb + 2 j. 


VIII. SUMMARY AND CONCLUSIONS 

The gauged B — L model is well-motivated and phenomenologically rich. Three heavy 
neutrinos, required for anomaly cancellation, participate in the seesaw mechanism and 

4 Same sign dileptonic signature can also arise from a Higgs triplet Type-II seesaw scenario [129] 
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generate viable light neutrino masses and mixings. In addition to the heavy neutrinos, 
this model has an extra Z' gauge boson and an additional SM-singlet Higgs state. In this 
work, we studied in detail the discovery prospect of the gauged B — L model through 
the Higgs window. We considered the B — L breaking scale to be few TeV, for which the 
Z' and the heavy neutrinos are naturally of TeV scale, and we considered heavy Higgs 
masses starting from 250 GeV up to 900 GeV. 

The SM-singlet Higgs state mixes with the SM Higgs with a mixing angle 9, that 
is constrained by direct and indirect searches, such as, vacuum stability, electroweak 
precision data and My/ mass measurement. The precise determination of the SM Higgs 
masses and the compatibility of its couplings with the SM prediction bounds the mixing 
angle to be sin 2 6 < 0.3. However, the most stringent limit on the mixing angle arises 
from the My/ mass measurement. In section IV, we reviewed the different constraints 
and for our subsequent analysis we consider the mixing angle sin 6 = 0.2, compatible with 
My/ mass measurement. 

During future LHC runs, Higgs coupling measurements can provide strong indirect 
constraints on the mixing angle between the SM-like Higgs boson and the heavy neutral 
Higgs, independent of the mass of the heavy state. The predicted bound on sin 6 from the 
14 TeV run of LHC with f Cdt = 300 fb _1 is sin 6 ~ 0.36 from the H —» WW* coupling 
measurement. However, if m# 2 < 500 GeV direct searches will be more sensitive. 

The heavy Higgs of this model can be produced via gluon fusion, VBF, associated 
Higgs production, out of which the gluon fusion offers the highest cross-section. We 
considered pp —* H 2 folded with H 2 decaying into different SM states. For our parameters 
of interest, H 2 does not decay to any additional invisibles state. The produced H 2 decays 
predominantly to W + W~, H\H\ and ZZ final states. The decays of the gauge bosons 
lead to different final states, such as 4£ and ljj$ T . We studied the discovery prospect of 
a heavy Higgs H 2 in the 4£, 2£2j and the £jj$r channels at the LHC (with f Cdt = 100 
fb _1 ) and HL-LHC (J Cdt = 3000 fb -1 ), where we employed a boosted decision tree to 
separate signal from background. 

The channel with four leptons was found to be the cleanest. The signal and background 
cross-sections for these processes are as — 0.1 fb and ctb — 42 pb, respectively. Using the 
cuts on the i) invariant mass of 4£ and on the reconstructed Z bosons, ii) the px cuts on 
the momenta of four leptons, as well as, the reconstructed Z bosons, we found that for a 
mass Mh 2 < 500 GeV, the H 2 can be discovered with a significance of ~ 5 a at HL-LHC 
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with 3000 fb _1 . 

The ZZ —>■ 212 j final state has a larger cross section than the 4£ final state. Partic¬ 
ularly for heavy H 2 masses an increased cross section is important to extend the LHC’s 
reach. However, we find that for small mixing angles (sin# = 0.2) this channel has small 
S/B and sensitivity. 

pp —> H 2 —;> W + W~ —y £jj$x offers an even larger cross-section for the signal, i.e. 
a ~ (9(10) fb. Unfortunately, huge backgrounds Ob — (9(10 3 ) pb make this channel 
extremely difficult to deal with. We explored the different mass regions of the heavy 
Higgs Mh 2 and applied exclusive trigger and selection cuts. We found that in the most 
optimistic scenario, assuming zero systematic uncertainties, a heavy Higgs search for 
350 < M H2 < 700 GeV has a statistical significance of 2.24<x with 100 flo -1 and up to 5 a 
with 3000 fb _1 . Since entirely negligible systematic uncertainties are not realistic in this 
fairly complex final state with jets and missing-transverse energy, even using data-driven 
methods only, the 4 1 channel proves to be superior in discovering a heavy Higgs boson in 
the B — L model. 

The heavy Higgs H 2 decays to H\Hi with a branching ratio ~ 20% for Mh 2 > 400 GeV. 
Due to the small production cross section for H 2 —> H\H\ —» bbbb or H 2 —>■ H\H\ —> 6677 
searching for H 2 in these final states is very challenging. However, we showed that for 
M H2 up to 500 GeV, the production cross-section for pp —> H 1 H 1 is significantly enhanced 
with respect to the SM expectation. Hence, it might be interesting to look for this channel 
in this region of the parameter space in more details. 

In this analysis, we also briefly considered the Z' searches. We used the ATLAS search 
at 8 TeV for Z' in the dileptonic channel to recast Z' mass constraints for various values 
of the U{1)b~l coupling g' and the heavy neutrino masses, M N . 
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